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ABSTRACT 

Significant adjustments to the values of the cosmological parameters estimated 
from high-redshift Type la Supernovae data are reported, almost an order of magni- 
tude greater than previously found. They arise from the effects of weak gravitational 
lensing on observations of high-redshift sources. The lensing statistics used have been 
obtained from computations of the three-dimensional shear in a range of cosmological 
iV-body simulations, from which it is estimated that cosmologies with an underly- 
ing deceleration parameter q = —0.51 + 0.03/ — 0.24 may be interpreted as having 
qo = —0.55 (appropriate to the currently popular cosmology with density parameter 
Om = 0.3 and vacuum energy density parameter Q\ = 0.7). In addition, the standard 
deviation expected from weak lensing for the peak magnitudes of Type la Super- 
novas at redshifts of 1 is expected to be approximately 0.078 magnitudes, and 0.185 
magnitudes at redshift 2. This latter value is greater than the accepted intrinsic dis- 
persion of 0.17 magnitudes. Consequently, the effects of weak lensing in observations 
of high-redshift sources must be taken properly into account. 

Key words: galaxies: clusters: general — cosmology: miscellaneous — cosmology: 
gravitational lensing — large-scale structure of Universe — cosmology: observations 
— supernovae: general 



1 INTRODUCTION 
1.1 Background 

The weak gravitational lensing of light from distant sources 
by the large scale structure in the universe results in the 
appearance of shear and convergence in images. The appli- 
cation of 'ray-tracing' methods, and other programmes, to 
cosmological iV-body simulations has enabled weak lensing 
statistics to be recorded for imaginary sources at various 
redshifts. These statistics in general show significant differ- 
ences between different cosmological models, and, of course, 
differences for sources at different redshifts. 

From weak lensing magnifications computed from 'ray- 
tracing' in iV-body simulations, Wambsganss et al. (1997) 
were able to estimate the likely magnification biases for 
sources at redshifts z = 0.5 and z = 1, and suggested that 
these biases should be applied to the observed magnitudes 
of high-redshift Type la Supernovae. Their work was con- 
ducted in a cosmological simulation with density parameter 
Hm = 0.4 and vacuum energy density parameter J7a = 0.6, 
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for which the deceleration parameter is go = —0.4. They 
found for this cosmology that if the computed median de- 
magnification for z = 1 were applied to the observed Su- 
pernovas at this redshift it would shift their positions on 
the Hubble diagram indicative of a different cosmology with 
qo = -0.395. 

Riess et al. (1998) and Perlmutter et al. (1999) have 
independently used data sets of Type la Supernovae extend- 
ing to a redshift of 0.97 to estimate values of the cosmolog- 
ical parameters from the Hubble diagram. In discussion of 
their results in the light of the work by Wambsganss et al. 
(1997), Riess et al. (1998) have stated that the effects on 
values of the cosmological parameters should be negligible. 
Perlmutter et al. (1999) have assumed that the effects of 
magnification or demagnification average out, and that the 
most overdense (or high magnification) lines of sight should 
be rare for their set of 42 high-redshift Supernovas. 

However, Barber et al. (2000) have now applied a new 
algorithm developed by Couchman, Barber and Thomas 
(1999) for the three-dimensional shear to obtain weak lens- 
ing statistics in simulations with cosmological parameters 
similar to those estimated from the high-redshift Super- 
novae data. Their results show that the median demagnifica- 
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tion values are considerably less and the dispersions in the 
probability distributions for the magnification are greater 
than found by Wambsganss et al. (1997) for redshifts of 
1 and 0.5 (although from different cosmologies). This sug- 
gests that the work of Wambsganss et al. (f997) should be 
reviewed in the light of the new weak lensing data now avail- 
able from a suitable cosmology. 

In this paper I apply the new data from Barber et al. 
(2000) in an attempt to re-estimate the effects of weak lens- 
ing on the cosmological parameters, and find the effects to 
be far from negligible. 

1.2 Outline of paper 

A brief outline of this paper is as follows. 

In Subsection 1.3 work by others on adjustments to 
the cosmological parameters resulting from weak gravita- 
tional lensing is summarised. 

Section 2 summarises the method followed by Bar- 
ber et al. (2000) for obtaining the weak lensing statistics. 
Subsection 2.1 describes the cosmological iV-body simu- 
lations; Subsection 2.2 summarises the main features of 
the algorithm used for the determination of the shear in 
three-dimensions; Subsection 2.3 outlines the method for 
integrating the shear values throughout the linked simula- 
tion volumes, and describes how these values are then used 
to determine the probability distributions for the magnifi- 
cations of sources at high-redshift. 

Section 3 summarises the relevant weak lensing statis- 
tics which are used in the present work. 

In Section 4, 1 describe the form of the Hubble diagram 
and its dependence on the cosmological parameters, and also 
summarise the important work by Riess et al. (1998) and 
Perlmutter et al. (1999) in estimating the cosmological pa- 
rameters from high-redshift Type la Supernovae data. 

In Section 5, I discuss how the dispersions in the mag- 
nification distributions for different source redshifts may af- 
fect the observed distance moduli. For the specific case of 
the 57m = 0.3, 57a = 0.7 cosmology, I then consider by how 
much the Supernovae data are displaced on the Hubble dia- 
gram, and how this translates into effective changes to the 
perceived values of the cosmological parameters. 

Section 6 summarises the results obtained for the de- 
celeration parameter, and discusses the findings in the light 
of work by others. 

1.3 Other work 

Wambsganss et al. (1997) have used the same ray-tracing 
method as described by Wambsganss, Cen and Ostriker 
(1998) for cosmological simulations with 57m = 0.4, 57a = 
0.6, and normalisation cr 8 = 0.79. The magnification values 
above and below which 97 \% of all of their lines of sight fall 
were /ii ow = 0.951 and /^high = 1.101 for source redshifts of 
z — 1, and /iw = 0.978 and /ihigh = 1.034 for z — 0.5. The 
median values for the magnification were 0.983 at z = 1, and 
0.993 at z = 0.5. The authors claim that these values would 
give rise to observed values of qo = —0.395 + 0.125/ — 0.095 
for the z = 1 data, and q = -0.398 + 0.048/ - 0.077 for 
the z = 0.5 data, rather than the assumed value of —0.4. 
(The quoted errors arise from the magnification ranges de- 
scribed above in the asymmetrical distributions.) They also 



report that the lensing-induced dispersions in their critical 
density cosmology are three times larger, but this cosmol- 
ogy uses a normalisation of as = 1.05 which overproduces 
the present-day rich cluster abundances. For this cosmology, 
they report magnifications up to 100, and correspondingly 
highly dispersed distributions, very much larger than those 
of Barber et al. (2000) for z = 3. The high magnifications 
almost certainly derive from the low value of their (fixed) 
softening scale. 

Fluke, Webster and Mortlock (1999) and Fluke, Web- 
ster and Mortlock (2000) have used a 'ray-bundle' method in 
which a discrete bundle of light rays is traced. The method 
allows a direct comparison between the shape and size of 
the bundle at the observer and at the source plane, so 
that the magnification, ellipticity and rotation can be de- 
termined straightforwardly. The cosmological models inves- 
tigated were similar to those of Barber et al. (2000) . How- 
ever, the particles were considered as point masses, so that 
very high magnification values could be achieved in princi- 
ple; however, to alleviate this situation, the authors do not 
include bundles which pass within \[2 of the Einstein radius 
of any particle. Because of their use of point particles, they 
use the empty beam approximation which gives rise to mag- 
nification probability distributions with minimum magnifi- 
cations of 1 (and therefore mean and median values greater 
than, or equal to, one), high magnification tails, and broad 
dispersions in magnification. For z = 1, the distribution 
in the magnifications for the SCDM cosmology (57m = 1, 
57a = 0) is clearly broader than that for the low density 
(LCDM) cosmology (57 M = 0.3, 57a = 0.7), which in turn 
is broader than the open (OCDM) cosmology (57m = 0.3, 
57a = 0). This is also true of Barber et al.'s (2000) data, 
although Barber et al. (2000) find a completely different or- 
der for the cosmologies at higher source redshifts. Because 
of the high magnification tails in Fluke, Webster and Mort- 
lock's (2000) distributions, they discard those high magnifi- 
cation lines of sight which occur with low probability before 
defining the /ii ow and ^ihigh values. 

These magnification probability distributions are then 
used by Fluke and Webster (2000) to examine the effects of 
the weak lensing dispersion on measurements of go in both 
the empty beam (more appropriate for point particles) and 
full beam limits. They consider whether the resulting shifts 
in the Hubble diagram can be fitted with various cosmolog- 
ical models. With input data from the z = 1 magnification 
distribution for the flat LCDM simulations, they find that 
qo = —0.53 + 0.96/ — 0.02 from the empty beam distribu- 
tion. To convert their empty beam magnification values to 
filled beam values, Fluke and Webster (2000) use a sim- 
ple scaling relationship. However, since the filled beam ap- 
proach is inappropriate to their method, the resulting val- 
ues should be treated with caution. They find the revised 
value qo = —0.61 + 0.89/ — 0.02; i.e., a much larger depar- 
ture than for the empty beam approximation, but a similar 
dispersion. Again the dispersion can be partly explained in 
terms of their use of point particles, although they do at- 
tempt to modify their distributions to account for the high 
magnification tails. 

It should be mentioned that Wambsganss et al. (1997) 
have assumed that the observed Supernovas have magnifica- 
tions of unity, and then have assessed departures from the 
best fit cosmology, rather than assuming that the best fit 
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Cosmology Qm ^A T erg No. of Box side 

particles (h~ 1 Mpc) 

SCDM 1.0 0.0 0.50 0.64 128 3 100 

TCDM 1.0 0.0 0.25 0.64 128 3 100 

OCDM 0.3 0.0 0.25 1.06 86 3 100 

LCDM 0.3 0.7 0.25 1.22 86 3 100 



to form the continuous depiction of the universe back to the 
source redshift. 



2.2 The three-dimensional shear algorithm 

Couchman et al. (1999) describe in detail the algorithm for 
the three-dimensional shear which has been applied by Bar- 
ber et al. (2000) to the simulations for the different cosmolo- 
gies. The main features of the code are as follows. 

First, the algorithm uses a Fast Fourier Transform 
method in the Particle-Mesh part of the code, in which the 
density distribution is smoothed. The procedure makes use 
of the periodicity of the fundamental volume, so that the ef- 
fects of matter effectively stretching to infinity are included 
in the computations. Couchman et al. (1999) pointed par- 
ticularly to the importance of including the effects of matter 
beyond a single period orthogonal to the direction of the line 
of sight. 

Second, in order to evaluate the shear components cor- 
rectly, the algorithm works with the 'peculiar gravitational 
potential,' which describes departures in the potential from 
homogeneity, and which results in net zero mass for the uni- 
verse. 

Third, a key feature of the algorithm is the variable par- 
ticle softening. The feature enables particles in low density 
regions to have extended softening, so that nearby evalua- 
tion positions for the shear register a density rather than a 
complete absence of matter. By contrast, particles in highly 
clustered regions are assigned low softening values, and a 
selected minimum value is introduced which limits the pos- 
sibility of strong lensing behaviour. The minimum value of 
the softening in the work by Barber et al. (2000) was selected 
to be 10 -3 in box units, equivalent to 0.1/t -1 Mpc at z = 0. 
The variable softening feature enables a much more realistic 
depiction of the large-scale structure within a simulation to 
be made. 

The output from the algorithm is the six independent 
components of the three-dimensional shear evaluated at a 
large number of positions within each simulation cube. 

2.3 Method 

A regular rectangular grid of 100 x 100 directions through 
each box was established. (Since there were only small de- 
flections, and the point of interest was the statistics of out- 
put values, each light ray was considered to follow one of the 
lines defined by these directions through the boxes.) Each 
'ray' was then connected with the corresponding line of sight 
through subsequent boxes in order to obtain the required 
statistics of weak lensing. In the application of the code, 
the shear was evaluated at 1000 positions along each of the 
lines of sight, forming a regular grid of evaluation positions 
in each box. 

The evaluation of the three-dimensional shear within 
the volume of the boxes, used the appropriate angular di- 
ameter distances at every evaluation position, representing a 
distinct advantage over two-dimensional (planar) methods. 
Furthermore, the variable softening facility within the al- 
gorithm lead naturally to the assumption that the universe 
may be described in terms of the filled beam approach, and 
this was the approach adopted by Barber et al. (2000) . This 



Table 1. Parameters used in the generation of the four different 
cosmological simulations. 

cosmology comes from the most likely demagnification at 
the peak of the magnification distributions, which is the ap- 
proach I take in this paper. Although Fluke and Webster 
(2000) use the median magnification values, the values are 
always close to unity and greater than, or equal to, unity for 
the empty beam case. 



2 METHOD 

2.1 The Hydra A^-body simulations 

The cosmological iV-body simulations used in the study by 
Barber et al. (2000) were provided by the Hydra Consortium 
and produced using the 'Hydra' iV-body hydrodynamics 
code (Couchman, Thomas and Pearce, 1995). Results are 
reported on simulations from four different cosmologies, re- 
ferred to as the SCDM, TCDM, OCDM and LCDM cosmolo- 
gies. Each of the simulations uses a Cold Dark Matter-like 
spectrum, and the parameters used in the generation and 
specification of these simulations are listed in Table |lj Qm 
and J1a are the present-day values of the density parame- 
ter and the vacuum energy density parameter respectively. 
The power spectrum shape parameter, T, is set to 0.5 in the 
SCDM cosmology, but the empirical determination (Peacock 
and Dodds, 1994) of 0.25 for cluster scales has been used in 
the other cosmologies. In each case, the normalisation, as, 
on scales of 8/i _1 Mpc has been set to reproduce the number 
density of clusters (Viana and Liddle, 1996); h is the Hubble 
parameter expressed in units of 100km s _1 Mpc _1 . 

In the SCDM and TCDM cosmologies, the number of 
particles was 128 3 , leading to individual dark matter parti- 
cle masses of 1.29 x lO 11 /!" 1 solar masses. In the low density 
universes, the number of particles was 0.3 times the number 
in the critical density universes, leading to the same individ- 
ual particle masses. The simulation output times were cho- 
sen so that consecutive simulation volumes could be snugly 
abutted; the side-dimensions were 100/i~ Mpc in every case. 

Since each time-output in a given simulation run is gen- 
erated using the same initial conditions, particular struc- 
tures (although evolving) occur at the same locations in all 
the boxes, and are therefore repeated with the periodicity 
of the box. To avoid these correlations, each simulation box 
was arbitrarily translated, rotated (by multiples of 90°) and 
reflected about each coordinate axis, before linking together 

t http:/ /hydra. mcmaster.ca/hydra/index. html 
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is different from the assumptions of many other workers who 
frequently use point particles, or a limited form of fixed soft- 
ening, and use the empty beam approximation. The two ap- 
proaches give rise to quite different expectations and results, 
the most obvious being the following. First, strong lensing 
can occur with particles of small radii, leading to high mag- 
nification tails in the probability distributions. Second, mag- 
nification distributions in the empty beam approximation 
all have minimum magnifications of /i m i n = 1, whilst in the 
filled beam, /i m i n < 1; this may alter the dispersions and the 
median values in the two distributions, essential to an under- 
standing of the magnitudes expected for high-redshift Type 
la Supernova?. Finally, the mean values for the magnifica- 
tions can be calculated from the respective angular diameter 
distances in the different cosmologies, for the empty beam 
approximation; however, the mean values in the filled beam 
approximation are always 1. 

The six independent second derivatives of the peculiar 
gravitational potential were calculated by the code at each of 
the selected evaluation positions throughout each box. These 
were then integrated in small steps along each line, forming, 
essentially, a large number of deflection sites through each 
simulation box. The integrated values formed the input data 
to establish the elements of the two-dimensional Jacobian 
matrix on each of the lines of sight for each of the deflec- 
tion sites. The Jacobian matrix elements were then used 
together with the multiple lens-plane theory (summarised 
by Schneider, Ehlers and Falco, 1992) and the appropriate 
angular diameter distances to obtain values for the magnifi- 
cation, source ellipticity, shear and convergence from each of 
the 10,000 lines of sight throughout all the simulation boxes 
linked back to the required redshift. 

In the procedure to obtain the two-dimensional second 
derivatives of the effective lensing potentials, it was assumed 
that (a) the angular diameter distances varied linearly be- 
tween the evaluation step points (although they were eval- 
uated exactly at each step point), (b) the angular diame- 
ter distances varied continuously through the depth of each 
simulation box, as they would in the real universe, (c) the 
shearing of light was weak, so that 'rays' were considered to 
follow the straight lines of sight defined by the grid of eval- 
uation positions within each simulation box, and (d) the 
accumulating component values within the developing Ja- 
cobian matrices were computed using the assumption that 
the shear was weak. 
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0.884 
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0.934 
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Table 2. Various magnification statistics for the different cos- 
mologies, as described in the text. 



spectively. The distributions are all broader in the SCDM 
cosmology when compared with the TCDM cosmology, due 
to its more clumpy character. For source redshifts of 4 the 
OCDM and SCDM cosmologies have very similar distribu- 
tions even though the angular diameter distance multiply- 
ing factors are larger in the OCDM cosmology. The most 
significant feature for our purposes here is that for high 
source redshifts the magnification distributions are broad- 
est in the LCDM cosmology (and the maximum values of 
the magnification are greatest here), but for lower source 
redshifts the width of the distribution falls below that for 
the SCDM and OCDM cosmologies. In the present work, I 
shall be mostly concerned with the magnification values for 
sources at z = 1.0 in the LCDM cosmology. 



3 WEAK LENSING IN DIFFERENT 
COSMOLOGIES 

Barber et al. (2000) have reported significant ranges of mag- 
nifications in all the cosmologies for nominal source redshifts 
of 0.5, 1.0, 2.0, 3.0 and 4.0. From the magnification distri- 
butions they have computed the median values, /i poa k, the 
values, ^i ow and /^high, above and below which 97^% of all 
lines of sight fall, and also the rms deviations from unity for 
the magnifications. All the values mentioned are displayed 
in Table which is reproduced from Barber et al. (2000) . 

It is interesting to compare the distributions in the dif- 
ferent cosmologies. Figure jl| and (also reproduced from 
Barber et al., 2000) show the magnification distributions 
for all the cosmologies for source redshifts of 4 and 1, re- 



4 DETERMINATIONS OF THE COSMOLOGY 
FROM TYPE IA SUPERNOVA 

The significant ranges in the magnifications (dependent on 
the cosmology) which might apply to distant sources have 
been summarised in the previous section. In the absence 
of magnification (or demagnification) from the large-scale 
structure, it would be possible to determine the cosmological 
parameters, £Im and f^A, from the departures from linearity 
in the Hubble diagram, provided 'standard candle' sources 
together with good calibration were available for measure- 
ment at high redshift. This is precisely the route taken by 
a number of authors, most importantly Riess et al. (1998) 
and Perlmutter et al. (1999), both of whom have used high- 
redshift Type la Supernova? data of redshifts up to 0.97. 
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Figure 1. The magnification probability distributions for all the 
cosmologies, assuming 2 = 4. 
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Figure 2. The magnification probability distributions for all the 
cosmologies, assuming 2 = 1. 



Carroll, Press and Turner (1992) give the general ex- 
pression for the distance measure to redshift, z: 



H D L 



— -sinn{| f2 k 



,1/2 



c(l + z) l^l 1 / 2 

f [(l + z') 2 (l + n M z')-z'(2 + z')Q A ]- 1/2 dz'}, (1) 
Jo 

where sinn means the hyperbolic sine, sinh, if Q,k > 0, or 
sine if Qk < 0; Ho is the Hubble parameter, Dl is the lu- 
minosity distance, related to the angular diameter distance 
through Dl = (1 + z) 2 D, and c is the velocity of light; 
the curvature density parameter, tlk, is defined in terms of 
the matter density parameter, Om, and the vacuum energy 
density parameter, SIa: 

l = Q M +^A+^ k . (2) 
However, if Qk = (which I have assumed for this work), 



hqd l r 

c(l + *) Jo 



[(1 + z')\\ + n M z') - z'(2 + z')n A ]- 1/2 dz' 



(3) 



and this distance measure, /, can be evaluated straightfor- 
wardly by numerical integration for any cosmology and any 
redshift. At low redshifts the expression reduces to the well- 
known Hubble law, but departures from the linear Hubble 
relation become evident as z is increased, and then different 
cosmologies may be differentiated. 

Observationally, the distance modulus, m — M, (where 
m is the apparent magnitude, and M the 'corrected' absolute 
magnitude) is determined, and this is related to the distance 
measure through 



m — M = 25 + 51og 10 [3000(1 + z)I] - 51og 10 /i; 



(4) 



see, e.g., Peebles (1993). 

Consequently, the value of the distance modulus is fixed 
for any given redshift, cosmology and Hubble parameter, but 
if the source is magnified (or demagnified) by weak gravi- 
tational lensing effects, then it will appear to be closer (or 
further away), according to 



Am = — 2.51og 10 /i, 



(5) 



where /i is the magnification. Alternatively, the observa- 
tional evidence could be interpreted as pointing to different 
values of the cosmological parameters, Qj» and Qa, since 
the curve on the Hubble diagram will be slightly displaced. 

I now briefly summarise the work of Riess et al. (1998) 
and Perlmutter et al. (1999), who have been able to estimate 
the cosmological parameters from data-sets of high-redshift 
Type la Supernova. 

Riess et al. (1998) have analysed the spectral and pho- 
tometric observational data of a total of 50 Type la Super- 
novae with redshifts 0.01 < z < 0.97. They have made use of 
two different correction methods for the peak magnitudes, 
namely the Multicolour Light Curve Shape method (MLCS) 
and a template method. The MLCS method, due to Riess, 
Press and Kirshner (1996) is an extension of Riess, Press 
and Kirshner's (1995) empirical method, described above, 
and makes use of up to 4 colours in the Supernovas photom- 
etry to quantify the amount of reddening by interstellar ex- 
tinction; in addition, Riess et al. (1998) have improved the 
method by using a larger 'training set' of Supernovae and 
including a quadratic term in the algorithm describing the 
light curve shapes. The template method is that of Hamuy 
et al. (1995) combined with Phillips' (1993) light curve de- 
cline parameter, Ami5(B). The fitting of the data has been 
done using a \ 2 statistic for a wide range of the parame- 
ters, Ho, Qm and f^A, and including both the dispersions in 
the host galaxy redshifts (in terms of the distance moduli), 
and the intrinsic dispersions (after correction) in the Super- 
novae peak magnitudes. With the constraint Qm + ^a = 1, 
Riess et al. (1998) find that Qa > at the 7a and 9a con- 
fidence levels for the two correction methods, and formally 
quote, with this constraint, 



Q M = 0.28 ±0.10, Ov = 0.72 ±0.10 

using the MLCS method, and 

n M = 0.16 ±0.09, Q A = 0.84 ±0.09 



(6) 



(7) 
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using the template method. 

Perlmutter et al. (1999) have based their analysis on 
42 high-redshift Type la Supernova? having redshifts 0.18 < 
z < 0.83, and combined their results with 18 Type la Super- 
novae from the Calan/Tololo Supernova Survey having red- 
shifts less than 0.1. Their method for determining the peak 
magnitudes relies on an empirical relationship between the 
light-curve width and the luminosity for the various tem- 
plate light-curves used. First, the time axis of the template 
is dilated by the factor (1 + z) to allow for the cosmological 
lengthening of the timescale. Then the observed light-curve 
is fitted to the template by introducing a 'stretch factor,' s, 
on the time axis. Finally, the peak magnitude is corrected 
to 

m CO rr = m + A corr (s), (8) 

where the correction term, A corr (s) is a simple monotonic 
function of s. The correction terms (different for different 
wavebands) were determined by Perlmutter et al. (1997) 
by relationship to the decline Amis, and the assessment 
of the luminosity by Hamuy et al.'s (1995) method. They 
claim that the residual dispersion in the peak magnitudes 
for all the Supernova?, after applying the correction based 
on the light-curve width-luminosity relationship, is within 
about 0.17 magnitudes. They state, however, that it is not 
clear whether the dispersion is best modeled as a normal dis- 
tribution in terms of the flux, or a log-normal distribution 
(Gaussian in magnitude terms). This point is worth bear- 
ing in mind if Monte Carlo simulations are performed to 
sample both the intrinsic dispersion and the magnification 
distribution arising from weak lensing. The preferred fitting 
of the data by Perlmutter et al. (1999) was performed by 
minimising \ 2 using the magnitude residuals in the Hubble 
diagram, from which they find: 

0.8Q M -0.6O\ ^ -0.2 ±0.1. (9) 

With the constraint, Qm + f^A = 1, they then find: 

Q M = 0.28±Hs (^statistical)!^ (systematic). (10) 

It is worth noting, however, that the above results are 
far from conclusive. Goodwin et al. (1999) for example, have 
pointed to the consistency of the data with a zero cosmologi- 
cal constant, provided the Hubble parameter has a higher lo- 
cal value compared with the global value; Riess at al. (1999) 
have pointed to the possible evolution of Type la Super- 
nova? light-curve shapes arising from differences in the rise- 
times, which may affect the estimation of the peak magni- 
tudes, although this has been questioned by Aldering, Knop 
and Nugent (2000). 



5 THE EFFECTS OF WEAK LENSING ON 
DETERMINATIONS OF THE 
COSMOLOGICAL PARAMETERS 

Both groups of workers, i.e., Riess et al. (1998) and Perl- 
mutter et al. (1999), point to cosmologies which are close to 
the fiM = 0.3, Qa = 0.7 cosmological simulation Barber et 
al. (2000) analysed in terms of weak lensing. I shall there- 
fore discuss the dispersions in magnification and the impact 
on determinations of the deceleration parameter, qo, with 
regard to this assumed cosmology. 
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0.99 
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0.993 
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1.005 




Table 


3. The various 


magnifications measured from the distri- 


butions for Type la Supernovae at the redshifts quoted in 


column 


1. All the values have been determined for the LCDM cosmology. 
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0.99 


0.062 


0.038 


-0.048/+0.050 


-0.190 
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0.88 


0.052 


0.034 
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0.039 
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-0.073 
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0.41 


0.015 


0.011 
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-0.055 
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0.29 


0.008 


0.005 
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-0.032 


0.039 


0.21 
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0.002 


-0.004/+0.004 


-0.019 


0.024 


0.10 


0.001 


0.000 


-0.001/+0.001 


-0.005 
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Table 4. Differences in the distance modulus, when compared 
with unmagnified sources, for the appropriate computed magnifi- 
cations from the distributions for Type la Supernovas sources at 
the redshifts given in column 1. The values in the final column 
represent the total dispersion in the distance modulus for sources 
with (de) magnifications of ^i ow and fihigh- All the values have 
been determined for the LCDM cosmology. 

By considering Type la Supernova? sources at ten evenly 
spaced redshifts between z = and z — 1, the resulting 
distributions of the magnifications show the expected in- 
creasing dispersions, which could be interpreted as obser- 
vational dispersions in the peak magnitudes. In Table ||, I 
give the values computed for /ii ow , /i pca k, the rms devia- 
tion in magnification and ^ihigh, and, in Table ^, what these 
values mean in terms of differences in the distance modu- 
lus compared with unmagnified sources. The final column 
in Table ^ gives the total difference in distance modulus be- 
tween the values derived from fj,\ ovl and /ihigh • The values for 
the changes in the distance modulus for appropriate values 
of the (de) magnification are shown in Figure ^. The asym- 
metry in the magnification distributions is clearly shown in 
terms of the distance modulus changes. 

Riess et al. (1998) have recognised that weak lensing 
of high-redshift Type la Supernova? can alter the observed 
magnitudes, and quote the findings of Wambsganns et al. 
(1997) i.e., that the light will be on average demagnified by 
1.7% at z — 1 in a universe with a non-negligible cosmo- 
logical constant. They state that the size of this effect is 
negligible. Perlmutter et al. (1999) assume that the effects 
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Figure 3. Differences in the distance modulus for appropriate 
values of the magnifications, as computed from the distribution 
curves for Type la Supernova? sources at redshifts < z < 1 in 
the LCDM cosmology. The solid curves show the dispersion in dis- 
tance modulus arising from the rms fluctuations in magnification; 
the dotted line shoes the departure in distance modulus arising 
from /ipeak; the dashed line arises from /ti ow ; and the dotted and 
dashed line arises from /%igh- 

of magnification or demagnification will average out, and 
that the most over-dense lines of sight should be rare for 
their set of 42 high-redshift Supernova?. However, they note 
that the average (de) amplification bias from integration of 
the probability distributions is less than 1% for redshifts of 
2 < 1. 

My work has shown that, in the LCDM cosmology, a 
source at z — 1 will have a median deamplification of 3.4% 
(fipenk = 0.966), and 1.3% at z = 0.5. Whilst these demag- 
nification factors are still small, they can result in reposi- 
tioning of the Supernova? data on the Hubble diagram with 
significantly revised values of the cosmological parameters. 
In addition, there is a significant probability of observing 
highly magnified Supernova?; 97|% of all lines of sight dis- 
play a range of magnifications up to 1.191 at z = 1, and 
this range is equivalent to a dispersion in the magnitudes of 
0.252. The standard deviation in the asymmetrical distribu- 
tions at z = 1 is 0.078 magnitudes. This is to be compared 
with the accepted intrinsic dispersion of 0.17 magnitudes for 
Type la Supernova? reported by Hamuy et al. (1996) in a 
set of low-redshift Supernovae. For sources at a redshift of 2, 
however, Barber et al.'s (2000) data predict a standard de- 
viation of 0.185 magnitudes from weak lensing, in excess of 
the intrinsic dispersion. Clearly, the weak lensing dispersion 
will become increasingly important as Supernovae at greater 
redshifts are discovered. 

In the work by Wambsganss et al. (1997) in the cosmol- 
ogy with Qm = 0.4 and Qa = 0.6 they found /ii ow = 0.951 
and /ihigh = 1.101 for source redshifts of z — 1, and /ii ow = 
0.978 and /ihigh = 1.034 for z — 0.5. Their median magnifica- 
tion values were 0.983 for z = 1 and 0.993 for z — 0.5. These 
may be compared directly with the values shown in Table ^ 
for the magnifications in the Q = 0.3, Qa = 0.7 cosmology 
I have investigated, i.e., /ti ow = 0.944 and /ihigh = 1.191 for 
source redshifts of z = 1, and /xi ow = 0.981 and /ihigh = 1.070 



Figure 4. The distance modulus for z = 1 vs. !!m, with the 
constraint Qm + Qa = 1. 



for z = 0.5. The median values are 0.966 for z = 1 and 0.987 
for z = 0.5. 

Since the most likely effect of weak lensing is a slight 
demagnification of the source, it is probable that the origi- 
nating data used by Riess et al. (1998) and Perlmutter et al. 
(1999) may suffer from such a demagnification, and therefore 
should be re-positioned on the Hubble diagram at brighter 
positions (i.e., lower distance moduli). This would then point 
to a slightly different cosmology with a different deceleration 
parameter. I plot, in Figure kl the distance modulus values 
at z — 1, for various values of Qm, centred around = 0.3, 
and keeping Qm + Qa = 1. The plot has been generated by 
solving equation Q numerically. 

With no magnification, Qm = 0.3, Qa = 0.7, h = 0.65, 
and 2 = 1, the distance modulus is 44.2626. Now assume 
that this cosmology is the best fit to the observed data. 
Then, by making a redshift 1 Supernova brighter than ob- 
served (assuming it was observed at the peak of the magni- 
fication distribution, with /i pca k = 0.966) the distance mod- 
ulus would be reduced to 44.2250, equivalent to Qm = 0.33, 
and equivalent to changing qo from —0.55 to —0.51. For the 
rms deviations in the magnification, the distance modulus 
changes by +0.0476/ — 0.0500, equivalent to revised qo val- 
ues of —0.61 and —0.49. For the /ihigh and /ti ow values at 
z = 1, go changes to —0.75 and —0.48. 

The conclusion, on the evidence of these data, is that 
true underlying cosmologies with qo = —0.51 + 0.03/ — 0.24 
may be interpreted as having qo — —0.55, from the use of 
perfect standard candles (without intrinsic dispersion), aris- 
ing purely from the effects of weak lensing. This adjustment 
to the value of qo is almost an order of magnitude larger 
than that found by Wambsganss et al. (1997) based on a 
cosmology with Qm = 0.4, Qa = 0.6, and the dispersion in 
the values is slightly larger than their findings. 



6 SUMMARY AND CONCLUSIONS 

I have summarised the results for the magnification distri- 
butions for the different cosmologies for different source red- 
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shifts as obtained by Barber et al. (2000) in Table ^. At high 
redshift, the LCDM cosmology produces the highest mag- 
nifications, the broadest distribution curves, and the lowest 
peak values. For sources at z = 3.6 in the LCDM cosmol- 
ogy, 97§% of all lines of sight have magnification values up 
to 1.850. (The maximum magnifications, not quoted here, 
depended on the choice of the minimum softening in the 
code, although the overall distributions were very insensitive 
to the softening.) The rms fluctuations in the magnification 
(about the mean) were as much as 0.191 in this cosmology, 
for sources at z = 3.6. Even for sources at z = 0.5 there is a 
measurable range of magnifications in all the cosmologies. 

The immediate implication is the likely existence of a 
bias in observed magnitudes of distant objects, and a likely 
dispersion for standard candles, for example, Type la Su- 
pernova? at high redshift. In particular, the weak lensing 
dispersion from 97|% of the lines of sight expected in the 
peak magnitudes for Type la Supernovae at redshifts of 1 
may be as much as 0.252 magnitudes. The standard devi- 
ations in the distributions are 0.078 magnitudes for z = 1 
and 0.185 magnitudes for z = 2. These values are to be 
compared with the accepted dispersion of 0.17 magnitudes, 
so that we should expect to see an increasing dispersion for 
the peak magnitudes of Supernovae as they are discovered 
at higher redshifts. 

I have made use of the magnification statistics from the 
LCDM cosmology, to reanalyse the results for the cosmo- 
logical parameters determined from the high-redshift Type 
la Supernova? data. (Both Riess et al., 1998, and Perlmuter 
et al., 1999, point to cosmologies with parameters close to 
those of the LCDM model, i.e., Q M = 0.3, Q A = 0.7.) 

The dispersions in the magnification are monotonically 
increasing with redshift, and have rms fluctuations about the 
mean of 0.045 for z = 0.99, and 0.016 for z = 0.49. These 
have been translated into variations in distance modulus, 
which, in turn, suggest a different cosmological model for 
the data. On the assumption that Qm — 0.3, Qa = 0.7 and 
go = —0.55 is the cosmology which fits the observed data 
best, my results indicate that the true underlying cosmology 
could have qo = —0.51 based on the peak magnification 
values, with a spread of +0.03/— 0.24 based on the computed 
magnification dispersion. This would represent a significant 
adjustment to the cosmological parameters. 

Wambsganss et al. (1997) found a median demagnifica- 
tion value much closer to unity (1.7% below for z = 1) in 
their cosmology with Qm = 0.4, Qa = 0.6 than Barber et al. 
(2000) have found (3.4% below for z = 1) in the Q M = 0.3, 
Qa = 0.7 cosmology. Consequently, my adjustment to go 
is considerably greater. However, the magnification disper- 
sions give rise to similar dispersions in go. Wambsganss et 
al.'s (1997) calculated adjustment to go is approximately an 
order of magnitude smaller than the departure I would sug- 
gest. 

Fluke, Webster and Mortlock's (2000) estimate for the 
revised value of the deceleration parameter (go = 0.53 + 
0.96/ - 0.02 for the Q M = 0.3, Qa = 0.7 cosmology) in 
the empty beam approximation is only half the value of the 
departure I suggest. The value of their dispersion can not be 
easily compared with mine, because of their use of the empty 
beam approximation, which gives rise to a different form for 
the distribution curve with values of magnification which 
are always positive. Furthermore, their method to convert 



the empty beam magnification values (more appropriate to 
their method) to filled beam values was not rigorous, so that 
the resulting values should be treated with caution. 

Both Wambsganss et al. (1997) and Fluke, Webster and 
Mortlock (2000) have assumed that the observed Super- 
novae have magnifications of unity, so that the median de- 
magnifications when using the filled beam approach are as- 
sumed to reposition the Supernovae at dimmer magnitudes. 
Fluke, Webster and Mortlock (2000) have to reposition them 
at brighter magnitudes when using the empty beam approx- 
imation. I have adopted the opposite viewpoint. I have as- 
sumed that the Supernovae are observed at their median de- 
magnification values, and therefore have to be repositioned 
at brighter magnitudes (i.e., reduced distance moduli) to ob- 
tain the correct cosmological parameters. Consequently, the 
direction of my adjustment to go is opposite to those of the 
above authors for the filled beam approximation. 
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